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Aims. Very-high-energy (VHE; ^100 GeV) y-rays are expected from y-ray bursts (GRBs) in some scenarios. Exploring this photon energy regime 
is necessary for understanding the energetics and properties of GRBs. 

Methods. GRBs have been one of the prime targets for the H.E.S.S. experiment, which makes use of four Imaging Atmospheric Cherenkov 
Telescopes (IACTs) to detect VHE y-rays. Dedicated observations of 32 GRB positions were made in the years 2003-2007 and a search for VHE 
\ y-ray counterparts of these GRBs was made. Depending on the visibility and observing conditions, the observations mostly start minutes to hours 
• r- I . after the burst and typically last two hours. 

■ Results. Results from observations of 22 GRB positions are presented and evidence of a VHE signal was found neither in observations of any 
!_h ' individual GRBs, nor from stacking data from subsets of GRBs with higher expected VHE flux according to a model-independent ranking scheme. 
Upper limits for the VHE y-ray flux from the GRB positions were derived. For those GRBs with measured redshifts, differential upper limits at 
the energy threshold after correcting for absorption due to extra-galactic background light are also presented. 

Key words, gamma rays: bursts - gamma rays: observations 



1. Introduction 

Gamma-ray bursts (GRBs) are the most energetic events in the y-ray 
regime. Depending on their duration (e.g. T 90 ), GRBs are categorized 
into long GRBs (Tg p > 2 s) and short GR Bs (Tgo < 2 s). First de- 
tected in late 1960s (Klebesadel et all 1 1973h . GRBs remained mysteri- 
ous for three decades. Breakthroughs in understanding GRBs came only 
after the discovery of lon ger-wavelength afterglo ws with the launch 
of BeppoSAX in 1997 ( Ivan Paradiis et all l2000h . Multi-wavelength 
(MWL) observations have proved to be crucial in our understanding 
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of GRBs, and provide valuable information about their physical prop- 
erties. These MWL afterglow observations are generally explained by 
synchr otron emission from shocked electrons i n the relativistic fireball 
model iPiranl [l999: Zhang & Meszaros, 2004). A plateau phase is re- 
vealed in man y of the Swifi PiKF light curves, the origin of which is 
still not clear (Zhang et al., 2006). Observations of GRBs at energies 
>10 GeV may test some of the ideas that ha ve been suggested to ex- 
plain the X-ray observations iFan et afl,l2008l) . 

In the framework of the relativistic fireball model, photons with en- 
ergies up to ~10 TeV or higher are expected from t he GRB afterglow 
phase (Zhang & Meszaros, 2004; Fan & Piran, 2008). Possible leptonic 
radiation mechanisms include forward-shocked el ectrons up- scattering 
self-emitted synchr o tron p h otons (SSC proce sses; iDermer et all l200d : 
Zhang & Meszarol 1200 ll : iFan et all l2008h or photons from other 
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shocked regions (Wans et all 1200 If) . Physical parameters, such as the 
ambient density of the surrounding material (n), magnetic field equipar- 
tition fraction (e s ), and bulk Lorentz factor (r bu i k ) of the outflow, may 
be constrained by obse rvations at these energies ( Wang~et"aTl 1200 ll ; 
iPe'er & Waxman.120051) . 

A possible additional contribution to VHE emission relates to the 
X-ray flare phenomenon. X-ray flares are fou nd in more than 50% o f 
the Swift GRBs during the afterglow phase 1 Chinca rini et a A 120071) . 
The energy fluence of some of them (e.g. GRB 050502B) is com- 
parable to that of the prompt emission. Mo st of them are clustere d 
at ~10 2 -10 3 s after the GRB (see Figure 2 in IChincarini et all [2007). 
while late X-ray flares (>10 4 s) are also observed; when these happen 
they can cause an increase in the X-ray flux of an order of magni - 
tude or more over the power-law temporal decay dCurran et all [2008). 
The cause of X-ray flares is still a subject of debate, but correspond- 
ing V HE y-ray flares from inverse-Compton ( IC) processes ar e pre- 
dicted dWang et all [200g ; iGalli & Pirol 120071 iFanetall [2008b . The 
accompanying external-Compton flare may be weak if the flare orig- 
inated b ehind the extern al shock, e.g. from prolonged central engine 
activity dFan et all 120081) . However, in the external shock model, the 
expected SSC flare is very strong at GeV energies and can be read- 
ily de tected using a VHE instrument with an energy threshold of ~ 100 
GeV dGalli & Piroll2008l) . such as the H.E.S.S. array, for a typical GRB 
at z~ 1 . Therefore, VHE y-ray data taken during an X-ray flare may help 
for distinguishing the internal/external shock origin of the X-ray flares, 
and may be used as a d iagnostic tool f or the late ce ntral engine activity. 

IWaxman & Bahcalll d2000l) and iMurase et all d2008l) suggest that 
GRBs may be sources of ultra-high-energy cosmic rays (UHECRs). In 
this case, 7r-decays from proton-y interaction may generate VHE emis- 
sion. The VHE y-ray emission produced from such a hadronic compo- 
nent is generall y expected to decay more s l owly than the le ptonic sub- 
MeV radiation (Bottcher & Dermer, 1998). Dermer (2007) suggests a 
combined leptonic/hadronic scenario to explain the rapidly-decaying 
phase and plateau phase seen in many of the SwiftPiKF light curves. 
This model can be tested with VHE observations taken minutes to hours 
after the burst. 

Mo st searches for VHE y-ray s from GRBs have o btained negative 
results dConnaughton et all 119971 ; I Atkins et all 120051) . There may be 
indications of excess phot on events from some observations, but these 
results are not conc l usive {A menomor i et allll99 6; Padilla et all 1 1991 : 
Atkin s et"ail l2OO0l : iPoirier et all 120031) . Curr ently, the most sen sitive 
detectors in the VHE y-ray regime are IACTs. lHoran et all 12007) pre- 
sented upper limits from 7 GRBs observed with the Whipple Telescope 
during the pre-Swift era. Upper limits for 9 GRBs with redshifts that 
were either unknown or >3.5 were also reported by the MAGIC col- 
laboration ( Alb ert et all 120071) . In general, these limits do not violate 
a power-law extrapolation of the keV spectra obtained with satellite- 
based instruments. However, most GRBs are now believed to originate 
at co smological distan ces, therefore absorption of VHE y-rays by the 
EBL (Nikishov, 1962) must be considered when interpreting these lim- 
its. 

In this paper, observations of 22 y-ray bursts made with H.E.S.S. 
during the years 2003-2007 are reported. They represent the largest 
sample of GRB afterglow observations made by an IACT array and 
result in the most stringent upper limits obtained in the VHE band. 
The prompt phase of GRB 060602B was observed serendipitously with 
H.E.S.S. The re sults of observations be fore, during, and after this burst 
are presented in Aharonian et al. (200§). 



2. The H.E.S.S. experiment and GRB observation 
strategy 

The H.E.S.S. arrajQ is a system of four 13m-diameter IACTs lo- 
cated at 1 800 m above sea level in the Khomas Highland of Namibia 
(23° 16' 18" S, 16°30'00" E). Each of the four telescopes is located at a 
corner of a square with a side length of 120 m. This configuration was 
optimized for maximum sensitivity to ~ 100 GeV photons. The effective 
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collection area increases from ~10 3 m 2 at 100 GeV to more than 10 nr 
at 1 TeV for observations at a zenith angle (Z.A.) of 20°. The system 
has a point source sensitivity above 1 00 GeV of ~ 1 .4x 1 0~ 1 1 erg cirr 2 s~ 1 
(3.5% of the flux from the Crab nebula) for a 5<r detection in a 2 h ob- 
servation. Each H.E.S.S. camera consists of 960 photomultiplier tubes 
(PMTs), which in total provide a field of view (FoV) of ~5°. This rela- 
tively large FoV allows for the simultaneous determination of the back- 
ground events from off-source po sitions, so that no dedicated off run is 
needed ( Aharonian et al., 2006c). The slew rate of the array is ~100° 
per minute, enabling it to point to any sky position within ~2 minutes. 
The H.E.S.S. array is currently the only IACT array in the Southern 
Hemisphere used for an active GRB observing programmfl 

Th e trigger system of the H.E.S.S. array is described in Funk et al. 
d2004l) . The stereoscopic technique is used, i.e. a coincidence of at least 
two telescopes triggering within a window of (normally) 80 nanosec- 
onds is required. This largely rejects background events caused by local 
muons that trigger only a single telescope. 

The observations reported here were obtained over the period 
March 2003 to October 2007. The observations of two GRBs in 2003 
were made using two telescopes while the system was under construc- 
tion. Before 2003 July, each of the two telescopes took data separately. 
Stereo analysis was then performed on the data, which requires co- 
incidence of events to be determined offline using GPS time stamps. 
After the installation of the central trigger system in 2003 July, the 
stereo multiplicity requirement was determined on-line. All observa- 
tions since 2004 have m ade use of the completed four-telescope array 
and the stereo technique (Aharonian et al., 2006c). 

Most of the data were taken in 28 minute runs using wobble mode, 
i.e. the GRB position is placed at an offset, 9 ff xt , of ±(f.5 or CP.7 (in R.A. 
and Deck) relative to the centre of the camera FoV during observations. 

Onboard GRB triggers distributed by the Swift satellite, as well 
as triggers from INTEGRAL and HETE-II confirmed by ground-based 
analysis, are followed by H.E.S.S. observations. Upon the reception 
of a GGS0 notice fr om one of these satellites (with appropriate indi- 
cationj^ that the source is a genuine GRB), the burst position is ob- 
served if Z.A. 545° (to ensure a reasonably low-energy threshold) dur- 
ing H.E.S.S. dark tim^J An automated program is running on site to 
keep the shift crew alerted of any new detected GRBs in real time. 
Depending on the observational constraints and the measured redshifts 
of the GRBs reported through GCN circular^ observations of the burst 
positions are started up to ~24 hours after the burst time, typically with 
an exposure time of a; 120 minutes in wobble mode. The remarkably 
nearby, bright GRB 030329 was an exceptional case. It was not ob- 
served until 11.5 days after the burst because of poor weather, which 
prohibited observation any earlier. 

3. The GRB observations 

Thirty-two GRBs were observed with H.E.S.S. during the period from 
March 2003 to October 2007. After applying a set of data-quality cri- 
teria that rejects observation runs with non-optimal weather conditions 
and hardware status, 22 GRB observations were selected for analysis 
and are described in this section. 

3.1. Properties of the GRBs 

For each burst, the observational properties as obtained from the trigger- 
ing satellite are shown in TableQ] These include trigger number, energy 
band, fluence in that energy band, and the duration of the burst (T90). 
Whenever there were follow-up observations in the X-ray, optical, or 

2 http://www.lsw.uni-heidelberg.de/projects/hess/HESS/grbs.phtml 

3 The Gamma ray bursts Coordinates Network, 
http://gcn.gsfc.nasa.gov/ 

4 which include, e.g., burst position incompatible with known 
sources, and a high signal-to-noise ratio of the burst 

5 H.E.S.S. observations are taken in darkness and when the moon is 
below the horizon. The fraction of H.E.S.S. dark time is about 0.2 

6 http://gcn.gsfc.nasa.gov/gcn3_archive.html 
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Table 1. Properties of GRBs observed with H.E.S.S. from March 2003 to October 2007. 



GRB 


Satellite 


Trigger 
number 


R.A.^ 


Decl. fl 


Error" 

(") 


Energy banc 
(keV) 


Fluence^ 
(10 erg cm 


1 90 

) (s) 


X^O'R 1 


: Rank^ 


071003 


Swift 


292934 


20 h 07 m 24 s .25 


+ 10°56'48'.'8 


5.7 


15-150 


830 


-150 


V 


vv 


1.604^ 


5 


070808 


Swift 


287260 


00°27 m 03.36 


+01 10 34. 8 


1.9 


15-150 


120 


-32 


V 


V • 


9 


070724A 


Swift 


285948 


01 h 51 m 13*96 


-18°35'40'.'l 


2.2 


15-150 


3 


~0.4 


V 


X X 


0.457£ 


21 


07072 IB 


Swift 


285654 


02 12 32.95 


-02 11 40. 6 


0.9 


15-150 


360 


-340 


V 


Vx 


3.626* 


10 


070721 A 


Swift 


285653 


00 12 39.24 


-28 22 Ou. 6 


2.3 


15-150 


7.1 


3.868 


V 


V • 


20 


070621 


Swift 


282808 


21"35 10U4 


-24 49 03. 1 


2 


15-150 


430 


33 


V 


X . 




1 


0706 12B 


Swift 


282073 


17" 26 m 54 s .4 


-08 45 08. 7 


4.7 


15-150 


168 


13.5 


V 


X . 




15 


070429A 


Swift 


277571 


1 _r f\\w a r> 

19 50 m 48.8 


-32 24 17. 9 


2.4 


15-150 


91 


163.3 


V 


V • 




3 


0704 19B 


Swift 


276212 


21 h 02 m 49!57 


-31°15'49'.'7 


3.5 


15-150 


736 


236.4 


V 


V ■ 




7 


070209 


Swift 


259803 


03 h 04 m 50 s 


-47°22'30" 


168 


15-150 


2.2 


0.09 


X 


X . 


0.314?^, 


22 


061110A 


Swift 


238108 


22 h 25 m 09!9 


-02 o 15'3077 


3.7 


15-150 


106 


40.7 


V 


V ■ 


0.758i~ 


11 


060526 


Swift 


211957 


15 h 31 m 18!4 


+00°17'11'.'0 


6.8 


15-150 


126 


298.2 


V 


V • 


3.2P~ 


8 


060505 


Swift 


208654 


22 h 07 m 04t50 


-27°49'57'.'8 


4.7 


15-150 


94.4 


-4 


V 


V • 


0.0889^, 


18 


060403 


Swift 


203755 


18 h 49 ra 21!80 


+08°19'45'.'3 


5.5 


15-150 


135 


30.1 


V 


X . 




16 


UjUoUl 


oWljt 


i4o jzz 


13 JO 


-Zl Jj 41 


i 


1 ^ i 

1 J— 1 JU 


j i 


1 Q A 


V 


Vx 


l.56j_ 


2 


050726 


Swift 


147788 


13 h 20 m lZ30 


-32°03'5078 


6 


15-150 


194 


49.9 


V 


V • 


13 


050509C 


HETE-II 


H3751 


12 h 52 m 53 s .94 


-44°50'04'.'l 


1 


2-30 


60 


25 


V 


vv 




19 


050209 


HETE-II 


U11568 


08 h 26 m 


+ 19°41' 


420 


30^100 


200 


46 




X . 




14 


041211B2 


HETE-II 


H3622 


06 h 43 m 12 s 


+20°23'42" 


80 


30^00 


1000 


>100 




X . 




4 


041006 " 


HETE-II 


H3570 


00 h 54 m 50 s :23 


+01°14'04'.'9 


0.1 


30^100 


713 


-20 


V 


vv 


0.716^ 


6 


030821 


HETE-II 


H2814 


21 h 42 m 


-44° 52 


O 


30^100 


280 


23 








17 


030329 


HETE-II 


H2652 


10 h 44 ra 49!96 


+21°31'17'.'44 


10~ 


30-400 


10760 


33 


V 


vv 


0.1687" 


12 



" R.A., Decl., and the positional errors (90% containment) were taken from GCN Reports (http://gcn.gsfc.nasa.gov/report_archive.html I for 
GRB 061 1 10A - GRB 071003 and GCN Circulars otherwise. 

* Fluenc e and T go data for GRB 050726 - GRB 0706 12B were taken from lSakamoto et al.1 i2008h except that T 90 of GRB 060505 was 
taken from lPalmer et al.l j2006h. Fluence and T 90 data of GRB 030329 and GRB 030821 were taken from lSakamoto et alj d2005h . and those of 
GRB 041006 from lShirasaki etall ( 120081) . Other data were taken from GCN Circulars and HETE pages (http://space.mit.edu/HETE/Bursts I. 

c X: X-ray, O: optical, R: radio; "V" indicates the detection of a counterpart, "x" a null detection, and "." that no measurement was reported 
in the corresponding energy range, from http://grad40.as.utexas.edu/grblog.php 
d The relative expected VHE flux for each GRB is ranked according to the empirical scheme described in Sect. 13.31 
' iPerlev et al.1 j2008l) 
/ ICucchiara et alJj2007l) 

s iMalesani et aD i2007h 

* Redshift of a candidate host galaxy Berg er & Foxl {2007). 
' iFvnbo et al.li2007h 

j IBerg er & Gladdersl d2006h 

* lOfek et alj j2006h 

' Redshift according to de Pasquale et al. (2007), based on afterglow modelling 

"' A lthough this burst was referred to as GRB 041211 in various GCN Circulars, the proper name GRB 04121 IB (e.g., in Pelangeon et al., 
2006) should be used to distinguish it from another burst, GRB 041211 A (=H3621) which occurred earlier on the same day (Pelangeon, A., 
private communication). 

" ISoderberg etall J2006I) 

" The position error of this burst is large, sec Fig. [3] 
" IStaneketalJj2003l) 



radio bands, whether a detection has occurred (denoted by a tick V) or 
not (denoted by a cross x) is also shown. If no observation at a given 
wavelength was reported, a dot (.) is shown. The reported redshifts (z) 
of 10 GRBs are also presented, of which 6 are lower than one. Two ob- 
served bursts, GRB 070209 and GRB 070724A, are short GRBs while 
the rest are long GRBs. The po pulation of short GRBs has a redshift 
distrib ution (IBerger et all 12 007) significantly less than that of the long 
GRBs ( Jakobsson et al., 2006). Therefore, on average they are likely to 
suffer from a lower level of EBL absorption. 



X-ray flares were detected from three of the GRBs in the H.E.S.S. 
sample. They occurred at 273s after the burst for GRB 050726, 284s 
for G RB 050801, and 2.6 x 10 5 s for GRB 070429A dCurran et all 
2008). Unfortunately, the flares occurred outside the time windows of 
the H.E.S.S. observations. 



3.2. H.E.S.S. observations 

For each burst, the start time, T start , of the H.E.S.S. observations after the 
burst is shown in Table[2] Since an observing strategy to start observing 
the burst position up to -24 hours after the burst time is applied, the 
mean T start is of the order of 10 hours. The (good-quality) exposure time 
of the observations using jV te ] telescopes for each burst is included. The 
mean ZA. of the observations is also presented. 

3.3. The ranking scheme 

As mentioned in the introduction, there is no lack of models predicting 
VHE emission from GRBs. However, the evolution of the possible VHE 
y-ray emission with time is model-dependent. To give an empirical, 
model-independent estimate of the relative expected VHE flux of each 
GRB (which also depends on r sta rtX it is assumed that: (1) the relative 
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VHE signal scales as the energy released in the prompt emission, taken 
as a typical energy measure of a GRB. Hence Fvhe k Fi5-i50kev where 
^i5-i50kev is the fiuence in the Swift/BAT band. For bursts not triggered 
by BAT, the measured fiuence is extrapolated into this energy band; (2) 
the possible VHE signal fades as time goes on, as observed in longer 
wavelength (e.g. X-ray) data. In particular, the VHE flux follows the 
average decay of the X-ray flux and therefore F V he 00 ^"ls-isokev x f 13 
where / denotes the time after the burst a nd 1.3 is the average X-ray 
afterglow late-time power-law decay index (Nousek et al., 2006). Since 
in most cases the exposure time of the observations is much shorter 
than r start (the start time of the corresponding H.E.S.S. observations af- 
ter the trigger), the expected flux at T stmt can be used as a measure of 
the strength of the VHE signal, and therefore of the relative possibility 
of detecting a VHE signal from that GRB. By setting t to r start , we have 

F V HE x ^15-150keV X T st ^ (1) 

The rank of each GRB according to equation [T] is shown in the last 
column in Table [T] Note that redshift information (available for only a 
few GRBs), and thus the corresponding EBL absorption, is not taken 
into account in the ranking scheme. 

4. Data Analysis 

Calibration of data, event reconstruction and rejection of the cosmic- 
ray backgro und (i.e. y-ray ev e nt sele ction criteria) were performed as 
described in Aharon ian et alj d2006cl) . which employs the techniques 
described bv lffillasl dl996l) . 

Gamma-like events were then taken from a circular region (on- 
source) of radius f? cut centred at the burst position given in Table[T] The 
background w as estimated using th e reflected-region background model 
as described in lBerge et al.1 d2007h . in which the number of background 
events in the on-source region (Nog) is estimated from Ji lcg j on off-source 
regions located at the same 8 ff xt as the on-source region during the 
same observation. The number of y-like events is given by N m - aN g 
where N on is the total number of events detected in the on-source region 
and a = 1 ln mgion the normalization factor. 

Independent analy ses of variou s GRBs using different methods and 
background estimates (Berge et al., 2007) yielded consistent results. 



4.1. Analysis technique 

Two sets of analysis cuts were applied to search for a VHE y-ray signal 
from observational d ata taken with thr ee or four telescopes. These are 
'standard' cuts I Ahar onian et al.Lr2006ch and 'soft' cutj^|(the latter have 
lower energy thresholds, as described in Aharo nian et all [2006a). For 
standard (soft) cuts, cut = 0.11° (9 mt = 0.14°). While standard cuts 
are optimized for a source with a power-law spectrum of photon index 
T = 2.6, soft cuts are optimized for a source with a steep spectrum (r = 
5.0), and have better sensitivity at lower energies. Since EBL absorption 
is less severe for lower energy photons, the soft-cut analysis is useful 
in searching for VHE y-rays from GRBs which are at cosmological 
dista nces. For example, the photon indices of two blazars PKS 2005- 
489 dAharonian etail|2005l) and PG 1553+1 13 dAharonian etaLll2008h 
were measured to be T i 4. 

An exception to this analysis scheme is GRB 030329. As the central 
trigger system had yet to be installed when this observation was made, 
a slightly different analysis technique was used. The description of the 
image a nd analysis cuts used f or the data from GRB 030329 can be 
found in lAharonian et al.1 d2005h . For GRB 030821, only the standard- 
cut analysis (for two-telescope data) was performed (see Sect. l5.4t . 

The positional error circle of most GRBs, with the exceptions of 
GRB 030821, GRB 050209, and GRB 070209, is small compared to 
the H.E.S.S. point spread function (PSF). The 68% y-ray containment 
radius, 6>68, of the H.E.S.S. PSF can be as small as ~3', depending on 
the Z.A. and # ffset of the observations, and the analysis cuts applied. The 
68% containment radius, 9 6S , of the observations of GRB 050209 and 
GRB 070209 is about 9' using standard-cut analysifl slightly larger 



than the corresponding error circles. Therefore, point-source analyses 
were performed for all GRBs except GRB 030821, the error box of 
which is much bigger than the H.E.S.S. PSF (see Sect. |5.4| for its treat- 
ment). 



4.2. Energy threshold 

The energy threshold, Efo, is conventionally defined as the peak in the 
differential y-ra y rate versus energy cur ve of a fictitious source with 
photon index T (Konopelk o et all! 1999ft . This curve is a convolution 
of the effective area with the expected energy spectrum of the source 
as seen on Earth. Such energy thresholds, obtained by the standard-cut 
analysis and the soft-cut analysis for each GRB observation, are shown 
in Table[2] assuming T = 2.6. The energy threshold depends on the Z.A. 
of the observations and the analysis used. The larger the Z.A., the higher 
is the energy threshold. Moreover, soft-cut analysis gives a lower value 
of E th than that of standard-cut analysis. Note that y-ray photons with 
energies below £ t h can be detected by the telescopes. 



4.3. Optical efficiency of the instrument 

The data presented were also corrected for the long-term changes in the 
optical efficiency of the instrument. The optical efficiency has decreased 
over a period of a few years. This has changed the effective area and en- 
ergy threshold of the instrument. Specifically, the energy threshold has 
increased with time. Using images of local muons in t he FoV, this effect 
in the calculation of flux upper limits is corrected (c.f. Aharonian et al., 
20061). 



5. Results 

No evidence of a significant excess of VHE y-ray events from any of 
the GRB positions given in Table Q] during the period covered by the 
H.E.S.S. observations was found. The number of on-source (N ou ) and 
off-source events (N ofl ), normalization factor (a), excess, and statis- 
tical significance of the excess in standard deviations (cr) are given 
for each of the 21 GRBs in Tabled The results for GRB 030821 are 
given in Sect. 15.4-1 Figure [T] shows the distribution of the significance 
obtained from the soft-cut analysis of the observations of each of the 
21 GRBs. A Gaussian distribution with mean zero and standard devi- 
ation one, which is expected in the case of no detection, is shown for 
comparison. The distribution of the statistical significance is consistent 
with this Gaussian distribution. Thus no significant signal was found 
from any of the individual GRBs. A search for serendipitous source 
discoveries in the H.E.S.S. FoV during observations of the GRBs also 
resulted in no significant detection. The 99.9% confidence level (c.l.) 
flux upper limits (above E h) have been calculated using the method of 
iFeldman & CousinsI d 19981) for both standard cuts (assuming T = 2.6) 
and soft cuts (assuming T = 5), and are included in Table[2] The limits 
are as observed on Earth, i.e. the EBL absorption factor was not taken 
into account. The systematic error on a H.E.S.S. integral flux measure- 
ment is estimated to be ~20%, and it was not included in the calculation 
of the upper limits. 

For those GRBs with reported redshifts, the effect of the EBL on 
the H.E.S. S. limits can be estimate d. Using the EBL model P0.45 de- 
scribed in Ahar onian et all d2006bl) . differential upper limits (again as- 
suming T = 5) at the energy threshold were calculated from the integral 
upper limits obtained using soft-cut analysis. These upper limits, as well 
as those calculated without taking the EBL into account, are shown in 
Table[3] 



7 'Soft' cuts were called 'spectrum' cuts in lAharonian et al] d2006al) . 

8 #68 is larger using soft-cut analysis 



calculated by eq. (17) in I.i <fc \t;j dl983h 



Table 2. H.E.S.S. observations of GRBs from March 2003 to October 2007. 



Standard-cut analysis 



GRET 




Exposure 


Mel Z.A. 


A'on 


A^OFF 


a Excess Signi- 


Efh 


Flux ULs 


M>N 


TVoFF 


a 


3xcess Signi- 


£th 


Flux ULs 


X 2 l&.o.f. 


POr 2 ) 




(min) 


(min) 




C) 










ficance (GeV) 


(cm ^ s 


) 










ficance (GeV) 


(cm s 


) 






070621 


6.5 


234.6 


4 


16 


204 


2273 


0.091 


-2.6 


-0.18 


250 


2.8 x 10~ 


2 


731 


5903 


0.13 


-6.9 


-0.24 


190 


5.6 x 10" 


2 


19.2/28 


0.89 


050801 


15.0 


28.2 


4 


43 


13 


173 


0.091 


-2.7 


-0.68 


400 


3.2 x 10" 


12 


46 


442 


0.13 


-9.3 


-1.2 


310 


1.6 x 10" 


1 


0.168/3 


0.98 


070429A 


64 


28.2 


4 


23 


4 


78 


0.091 


-3.1 


-1.2 


290 


2.4 x 10" 


12 


20 


203 


0.13 


-5.4 


-1.0 


220 


1.0 x 10" 


I 


6.39/3 


0.094 


04121 IB" 


1 567.1 


14.2 


3 


64 


9 


87 


0.11 


-0.67 


-0.21 


1850 


6.8 x 10~ 


2 


27 


236 


0.17 


-12 


-1.9 


1360 


2.6 x 10" 


1 


| 14.6/14 


0.40 


{ 742.3 


112.3 


4 


44 


76 


1247 


0.063 


-1.9 


-0.21 


380 


3.7 x 10~ 


2 


317 


4353 


0.083 


-46 


-2.4 


280 


1.8 x 10" 


1 


071003^ 


623.3 


56.2 


4 


35 


16 


272 


0.10 


-1 1 


-2.2 


390 


1.0 x 10~ 


12 


97 


785 


0.14 


-15 


-1.4 


280 


1.4 x 10" 


1 


> 

\ 32.3/12 


0.0012 


691.1 


56.2 


3 


41 


25 


204 


0.10 


4.6 


0.93 


480 


5.6 x 10~ 


12 


79 


547 


0.14 


0.86 


0.091 


340 


1.5 x 10" 


1 


041006 


626.1 


81.9 


4 


27 


80 


770 


0.10 


3 


0.32 


200 


1.1 x 10~ 


11 


302 


1974 


0.14 


20 


1.1 


150 


6.8 x 10" 


1 


8.89/9 


0.45 


0704 19B 


907 


56.4 


4 


47 


28 


391 


0.091 


-7.5 


-1.3 


700 


2.4 x 10" 


2 


121 


1069 


0.13 


-13 


-1.0 


520 


7.5 x 10" 


2 


11.9/6 


0.064 


060526 


284.2 


112.8 


4 


25 


93 


1068 


0.10 


-13.8 


-1.3 


280 


2.9 x 10" 


2 


492 


3711 


0.14 


-38 


-1.6 


220 


9.2 x 10" 


2 


19.8/12 


0.072 


070808 


306.2 


112.8 


4 


34 


49 


659 


0.091 


-11 


-1.4 


310 


3.2 x 10" 


12 


209 


1733 


0.13 


-7.6 


-0.49 


260 


7.5 x 10" 


2 


15.8/12 


0.20 


07072 IB 


925.7 


103.8 


4 


40 


59 


984 


0.063 


-2.5 


-0.31 


440 


1.4 x 10" 


12 


237 


2676 


0.083 


14 


0.89 


320 


8.8 x 10" 


2 


15.5/11 


0.16 


061110A 


407.68 


112.8 


4 


25 


76 


838 


0.093 


-1.9 


-0.21 


280 


4.3 x 10" 


12 


314 


2671 


0.13 


-20 


-1.0 


200 


8.4 x 10" 


2 


4.66/11 


0.95 


030329^ 


16493.5 


28.0 


2 


60 


4 


26 


0.14 


0.27 


0.13 


1360 


2.6 x 10" 


12 


















5.93/3 


0.12 


050726~ 


772.7 


112.8 


4 


40 


107 


1031 


0.083 


21 


2.1 


320 


7.1 x 10" 


2 


333 


2619 


0.11 


42 


2.3 


260 


3.4 x 10" 


1 


14.7/12 


0.26 


050209 


1208.5 


168.6 


4 


48 


104 


1096 


0.11 


-18 


-1.6 


480 


4.4 x 10" 


2 


528 


4204 


0.14 


-73 


-2.8 


340 


1.5 x 10" 


1 


36.3/18 


0.0065 


0706 12B 


901.7 


112.8 


4 


18 


104 


1190 


0.091 


-4.2 


-0.39 


240 


4.1 x 10" 


12 


415 


3233 


0.13 


11 


0.51 


180 


1.5 x 10" 


1 


4.87/12 


0.96 


060403 


820.4 


52.8 


4 


39 


33 


252 


0.091 


10 


1.9 


440 


4.8 x 10" 


12 


128 


875 


0.13 


19 


1.6 


320 


1.3 x 10" 


1 


10.4/6 


0.11 


060505 


1163 


111 


4 


42 


99 


837 


0.091 


23 


2.4 


520 


5.6 x 10" 


2 


339 


2740 


0.13 


-3.5 


-0.18 


400 


3.9 x 10" 


2 


22.1/12 


0.036 


050509C 


1289 


28.2 


4 


22 


31 


344 


0.083 


2.3 


0.41 


200 


1.7 x 10" 


11 


112 


965 


0.11 


4.8 


0.43 


150 


1.5 x 10" 





0.301/3 


0.96 


070721 A 


893.5 


112.8 


4 


30 


90 


1436 


0.059 


5.5 


0.58 


320 


6.5 x 10" 


12 


280 


3837 


0.077 


-15 


-0.86 


260 


1.3 x 10" 


1 


6.78/12 


0.87 


070724A 


927.5 


84.6 


4 


23 


73 


720 


0.091 


7.5 


0.88 


260 


7.3 x 10" 


12 


246 


2042 


0.13 


-9.3 


-0.55 


200 


1.0 x 10" 


1 


14.3/9 


0.11 


070209 


926.7 


56.4 


4 


41 


37 


444 


0.091 


-3.4 


-0.51 


480 


2.3 x 10" 


2 


185 


1442 


0.13 


4.8 


0.33 


370 


1.1 x 10" 


1 


5.35/6 


0.50 



Soft-cut analysis 



Temporal analysis 



8 The GRBs are listed in the order of the ranking scheme described in Sect. 13. 31 GRB 030821 is not listed, the results of which are given in Sect. 
* Three- and four-telescope data are presented. 



A slightly different analysis technique was used, see Sect. 14. II Soft-cut analysis is not available for this observation. 
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Secondly, combining the significance of the results from three selected 
subsets extracted from the whole sample was performed. The a priori 
selection criteria were to choose those GRBs with a higher expected 
VHE flux or a lower level of EBL absorption. The following require- 
ments were used to select three subsets: 



-10 12 3 4 

Significance of individual GRBs 

Fig. 1. Distribution of the statistical significance {histogram) as 
derived from the observations of 20 GRBs using soft-cut anal- 
ysis. The mean is -0.4 and the standard deviation is 1 .4. Each 
entry corresponds to one GRB. The solid line is a Gaussian func- 
tion with mean zero and standard deviation unity. 



5.1. Stacking analysis 

Although no significant excess was found from any individual GRB, co- 
adding the excess events from the observations of a number of GRBs 
may reveal a signal that is too weak to be seen in the data from one GRB, 
provided that the PSFs of the H.E.S.S. observations are bigger than the 
error box of the GRB positions (which is the case, see Sect. l4.lt . Firstly, 
stacking of all GRBs (except GRB 030821, which has a high positional 
uncertainty) in the sample was performed. This yielded a total of -157 
excess events and a statistical significance of -1.98 using the soft-cut 
analysis. Use of standard cuts produced a similar result (see Table |4j. 



Table 3. Differential flux upper limits at the energy thresholds 
from the H.E.S.S. observations of GRBs with reported redshifts. 



GRB 


Redshift 


Ea, (GeV) 






r corrected 


060505 


0.0889 


400 


3.9x10" 


14 


5.8x10" 


14 


030329 


0.1687 


1360 


7.6x10" 


15 


9.7x10" 


14 


070209 


0.314 


370 


1.2x10" 


13 


8.7x10" 


13 


070724A 


0.457 


200 


2.1x10" 


13 


1.0x10" 


12 


041006 


0.716 


150 


1.8x10" 


12 


2.7x10" 


1 1 


061110A 


0.758 


200 


1.7x10" 


13 


1.7x10" 


11 


050801 


1.56 


310 


2.1x10" 


13 


b 




071003^. 


1.604 


280 


2.0x10" 


13 


T 




060526~ 


3.21 


220 


1.7x10" 


13 


b 




07072 IB 


3.626 


320 


1.1x10" 


13 


b 





.... ,„ „jT 2 s" 1 GeV" 1 . 

* The limits corrected for EBL absorption are >10 orders of magni- 
tude larger than that observed. 
Only 4-telescope data were used. 

Table 4. Combined significance of 3 subsets of GRBs selected 
based on the requirements listed in Sect. 15.11 





Number 


Soft-cut 


Standard-cut 




of GRBs 


analysis 


analysis 


Sample A 


10 


-2.13 


-1.81 


Sample B 


6 


-0.20 


1.45 


Sample C 


11 


-0.53 


0.48 


all GRBs 


21 


-1.98 


-0.18 



Sample A: the first 10 in the ranking described in Sect. 13.31 

Sample B: all GRBs with a measured redshift z < 1; 

Sample C: all GRBs with a soft-cut energy threshold lower than 

300 GeV and with either a measured redshift z < 1 or with an 

unknown redshift. 

The result is shown in Table|4] As can be seen, there is no significant 
evidence of emission in any of these subsets. 

5.2. Temporal analysis 

As possible VHE radiation from GRBs is expected to vary with time, 
a temporal analysis to search for deviation from zero excess in the ob- 
served data was performed. Soft-cut analysis was used for all GRBs (ex- 
cept GRB 030329) since this analysis has a lower energy threshold and 
a better acceptance of y-rays and cosmic rays and therefore increases 
the statistics. The y-like excess events were binned in 10-minute time 
intervals for each GRB data set and were compared to the assumption 
of no excess throughout the observed period. The^ 2 /d.o.f. value and the 
corresponding probability are shown in Table [2] for each GRB. Within 
the whole sample, the lowest probability that the hypothesis that the ex- 
cess was zero throughout the observation period is correct is 1.2 x 10" 3 
(for GRB 071003) and no significant deviation from zero within any 
of the GRB temporal data was found. Standard-cut analysis produced 
consistent results. 



5.3. GRB 070621: Observations of a GRB with the fastest 
reaction and the longest exposure time 



GRB 070621 is the highest-ranked GRB in the sample (Sect 133b. i.e. 
it has the highest relative expected VHE flux at the start time of the 
observations. The duration of the Swift burst was Tgo ~ 33s, thus 
clearly classifying the burst as a long GRB. The fluence in the 15- 
150 keV band was ~4.3xl0" 6 erg cm" 2 . The XRT light curve is rep- 
resented by an initial rapidly-decaying phase and a shallow phase, with 
the tr ansition happening aro und t + 380s where t denotes the trigger 
time dSbarufatti e t al., 2007). Despite extensive optical monitoring, no 
fading optical counterpart was found. The H.E.S.S. observations started 
at to + 420s and lasted for ~5 hours, largely coincident with the X-ray 
shallow phase. These observations were both the most prompt and the 
longest among those presented. Figure [2] shows the 99.9% H.E.S.S. en- 
ergy flux upper limits above 200 GeV ( using soft-cut analysis), together 
with the XRT results ( Evans et al 

1E007|). 

As seen, the limits for this pe- 
riod are at levels comparable to the X-ray energy flux during the same 
period. Unfortunately the lack of redshift information for this burst pre- 
vents further interpretation of the limits. 



5.4. GRB 030821: Observations of a GRB with a high 
positional uncertainty 

Some GRBs, such as GRB 030821, have a high uncertainty in position; 
with a relatively large camera FoV (~5°), the H.E.S.S. telescopes are 
able to cover the whole positional error box of such GRBs. 

Observations of GRB 030821 started 18 hours after the burst and 
lasted for a live-time of 55.5 minutes, with a mean Z.A. of 28°. The 
observations were taken when the array was under construction and 
only two telescopes were operating, resulting in an energy threshold 
of 260 GeV. The GRB has a relatively high uncertainty in position 
as de termined from IPN (the third Interplanetary Network) triangula- 
tion dHurlev et all 1200 3h . and its error box is bigger than the PSF of 
H.E.S.S. However, because of the relatively large FoV of the camera, 
the whole error box, and thus the possible GRB position, is within the 
H.E.S.S. FoV. The sky excess map overlaid with the error box is shown 
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6. Discussion 
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Fig. 2. The 99.9% confidence level energy flux upper limits (in 
red) at energies >200 GeV derived from H.E.S.S. observations 
at the position of GRB 070621. The ends of the horizontal lines 
indicate the start and end times of the observations from which 
the upper limits were derived. The XRT ener gy flux in the 0.3-1 
keV band is shown in black for comparison (Eva ns et aUl2007l) . 




21*50" 



2f40° 



21 n 35 m 
Right Ascension 



Fig. 3. The y-like excess events in the region of the GRB 03082 1 . 
The error box shows the position of the burst localized by IPN 
triangulation dHurley et al.L l2003). The colour (grey) scale is set 
such that the blue/red (black/grey) transition occurs at the ~ 1 .5<t 
significance level. The sky map was derived using two obser- 
vations pointing at two different positions (marked by crosses), 
resulting in a non-uniform distribution of events in the map. 



in Fig. [3] As can be seen, there is no significant excess at any position 
within the error box. The sky region with the largest number of peak 
excess events is located in the south-eastern part of the error box. Using 
a point-source analysis centred at this peak, a flux upper limit (above 
260 GeV) of ~1.7xl0~" cirT 2 s was derived. Since an upper limit 
derived for any location in the error box with fewer excess events is 
lower than this valuJ^T. it may be regarded as a conservative upper limit 
of the VHE flux associated with GRB 030821 during the period of the 
H.E.S.S. observations. 



The upper limits presented in this paper are among the most stringent 
ever derived from VHE y-ray observations of GRBs during the af- 
terglow period. In fact, the 99.9% confidence level limits (in energy 
flux) are at levels comparable to the X-ray energy flux as observed 
by SwiftPiKF during the same period (see, e.g. Fig. |2j. Unless most 
of the GRBs are located at high redshifts and thus their VHE flux is 
severely absorbed by the EBL (this possibility is discussed below), one 
expects detection of the predicted VHE component wi th energy flux 
levels comparable to those in X-rays in some scena rios (Dermer et al., 
2000 ; iWanget all 120011 , IZhang & Meszarosl. 1200 ll ; iPe'er & Waxmaii , 
booalFan et al.Ll2008h . 

On the other hand, the unknown redshifts of many of the GRBs 
in the sample (including GRB 070621, the highest-ranking, which is 
discussed in Sect. l3.3t complicate the physical interpretation of the data, 
because EBL absorption at VHE energies is severe for a GRB with z > 
1 . The mean and median redshift of the 10 GRBs with reported redshifts 
is 1.3 and 0.7, respectively. If the 12 GRBs without redshift have the 
same redshift distribution, one would expect ~40% of them (~5 GRBs) 
to have z < 0.5. In this case, the EBL absorption may not preclude the 
detection of the predicted VHE y-rays for the GRB sample presented 
hertO 

There is no reported X-ray flare during the H.E.S.S. observational 
time windows, therefore no conclusion on whether or not X-ray flares 
are accompanied by VHE flares, as well as the origin of X-ray flares, 
can be drawn. If UHECRs are generated in nearby GRB sources, as 
suggested by some authors, a detectable VHE flux is expected from 
nearby GRBs. Therefore, although the unknown redshifts of a signifi- 
cant fraction of GRBs in our sample (12 out of 22) and the uncertainty 
in the modeled VHE temporal evolution are surely in play, the results 
presented here do not indicate (but also not exclude) that GRBs are 
dominant sources of UHECRs. 



7. Outlook 

The data from our sample of 22 GRBs do not provide any evidence 
for a strong VHE y-ray component from GRBs during the afterglow 
phase. EBL absorption can explain the lack of detection in our sample. 
However, this does not exclude a population of GRBs that exhibit a 
strong VHE component. While the EGRET experiment did not detect 
MeV-GeV photons from most BATSE GRBs in its FoV, some strong 
bursts (e.g. GRB 940217) have proved to emit delayed em ission, ~1.5 
hours after the burst, at energies as high as ~20 GeV (Hurlev et al., 
1994). With Fermi's observations of GRBs having started in mid-2008, 
it is likely that our knowledge of the high-energy emission of GRBs will 
be improved in the near future. 

The future prospects for detection at VHE energies rely on the like- 
lihood of observing a GRB with low redshift (e.g. z < 0.5) early enough. 
In the cases where there is no detection, sensitive and early upper limits 
on the intrinsic VHE luminosity of these nearby GRBs will still im- 
prove our understanding of the radiation mechanisms of GRBs. The ex- 
istence of a distinct population of low-luminosity (LL) GRBs was sug- 
gested based on the high detection rate o f low-redshift LL GRBs such 
as GRB 980425 and GRB 060218 (e.g. ISoderberg et all 120041) . Due 
to their proximity, they are good targets for VHE observations. Since 
they are sub-energetic compared to other GRBs, they may be accompa- 
nied by a lower VHE luminosity. On the other hand, if most radiation 
are emitted at high energies, the detection probability would be much 
higher. 

iFranceschini et al] (12008) claimed a very small y-ray opacity due 
to EBL absorption. The optical depth is about a factor of three less tha n 
the one we used, depending on the energies ( Aharonian et al., 2006b). 
Therefore, on-going GRB observations with H.E.S.S., as well as other 
ground-based VHE detectors, are crucial to test this model. 



10 A larger excess implies a higher value of the upper limit, since the 
integrated exposure, that depends on ZA. and 9 a ff xt of the observations, 
is largely the same over the whole error box. 



11 The optical depth of EBL absorption for a ~ 100 G eV photon is ~3 
at z = 1 , according to the P0.45 model demonstrated in Aharonian et al. 
(20063). 
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8. Conclusions 

During 5 years of operation (2003-2007), 32 GRBs were observed dur- 
ing the afterglow phase using the H.E.S.S. experiment. Those 22 GRBs 
with high-quality data were analysed and the results presented in this 
paper. Depending on the visibility and observing conditions, the start 
time of the observations varied from minutes to hours after the burst. 

There is no evidence of VHE emission from any individual GRB 
during the period covered by the H.E.S.S. observations, nor from stack- 
ing analysis using the whole sample and a priori selected sub-sets of 
GRBs. Fine-binned temporal data revealed no short-term variability 
from any observation and no indication of VHE signal from any of these 
time bins was found. Upper limits of VHE y-ray flux during the obser- 
vations from the GRBs were derived. These 99.9% confidence level 
energy flux upper limits are at levels comparable to the contemporary 
X-ray energy flux. For those GRBs with reported redshifts, differential 
upper limits at the energy threshold after correcting for EBL absorption 
are presented. 

H.E.S.S. phase II will have an energy threshold of about 30 GeV 
With much less absorption by the EBL at such low energies, it is hoped 
that the H.E.S.S. experiment will enable the detection of VHE y-ray 
counterparts of GRBs. 
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